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Показана возможность использования решетчатых периодических функций при обработке узкополосных сигналов методами 
дискретного преобразования Фурье, мгновенной спектральной плотности и синхронного детектирования, обеспечивающая со- 
кращение количества обрабатываемых отсчетов. 


Задача определения спектрального состава: ам- 
плитуд, частот и фаз составляющих периодическо- 
го сигнала - встречается при анализе многочастот- 
ных колебаний в системах управления [ 1 ] , в задачах 
распознавания многочастотных сигналов в цифро- 
вых каналах связи [2], при разработке и эксплуата- 
ции систем диагностирования электротехнических 
и электромеханических устройств [3]. 

В настоящее время в качестве математического 
аппарата цифрового спектрального анализа широ- 
ко используют методы дискретного преобразова- 
ния Фурье, часто в виде быстрого преобразования 
Фурье, мгновенной спектральной плотности и 
синхронного детектирования. Однако реализация 
классических алгоритмов данных методов требует 
больших вычислительных затрат, что существенно 
усложняет техническую реализацию устройств и 
ограничивает их быстродействие. Стремление уме- 
ньшить объем вычислений стимулировало разви- 
тие новых алгоритмов и методов спектрального 
оценивания, направленных на устранение указан- 
ных недостатков. 

В настоящей работе рассматривается возмож- 
ность использования в качестве опорных, взамен 
дискретизируемых гармонических функций вида 
8Іп(ю/) и С08(ю/), решетчатых функций 
Аі„(4) = 8іп(2^4)=8іп(ю ( 4) и 4 к (4)=со8(2я^4)=со8(ю,4), 
принимающих на своем периоде строго фиксиро- 
ванное количество ( М) значений в определенные 
моменты времени і ь когда: 


"у ч 


2 тг 72. к 2 к . , ч 2к 

— , 1 17’ 2 77’-’ (М - 1) 17’ 

МММ м 


при различных значениях круговых частот опорно- 
го сигнала <у=2я/ где /- частота опорного сигна- 
ла (Гц), у - порядковый номер спектральной соста- 
вляющей. 

Рассмотрим особенности использования ре- 
шетчатых функций в спектральном анализе перио- 
дических сигналов методами: 

• мгновенной спектральной плотности (МСП) [1] ; 

• синхронного детектирования (СД) [5]; 

• дискретного преобразования Фурье (ДПФ) [4]. 


Мгновенная спектральная плотность 


Мгновенные спектральные плотности из [1] по- 
зволяют найти частоту той или иной составляющей 
сигнала по экстремуму 5(оо): 


Аф®) = |/(0 8Іп(й> І)с1ѵ, 

О 

Т 

5 2 (со) = |/(7)со8(й>Т)Л; >, 

О 

5(со) = ^5?(со) + 5; (со) , 


( 2 ) 


где Т - период сигнала, //) - анализируемый сиг- 
нал, со - круговая частота. 

Анализ ведется в определенном диапазоне ча- 
стот [с0 НЛЧ ,...,(О П р Е д]. 

При анализе дискретизированного сигнала, 
представленного мгновенными цифровыми отсче- 
тами, формулы (2) примут вид: 

/(«) = і /(7,.)5Іп(ю (.); 

1=0 
**и - 1 

/>(«,) = 1 /( 4 ) 008 ( 6 ^ 4 );^, ( 3 ) 

5(®,) = л /5 1 2 (®,)+ $1(0), ) , 

где с о ш <(о<(о П рщ , Ы и — количество отсчетов на ин- 
тервале дискретизации И и =Т и /Аі, Т и и А! - интер- 
вал и шаг дискретизации,//,) - мгновенное значе- 
ние анализируемого периодического сигнала, 
І - порядковый номер спектральной составляю- 
щей, 4 - дискретный момент времени. 

Зафиксировав частоты, на которых наблюдает- 
ся экстремум 5(со), можно провести оценку спек- 
трального портрета анализируемого периодическо- 
го сигнала. По значению 5(со) для каждой круговой 
частоты ю,=2я^ можно определить амплитуду соот- 
ветствующей спектральной составляющей [1]: 

А,у = №;)]№■ ( 4 ) 

Фазовый угол выявленной спектральной соста- 
вляющей рассчитывается по формуле [1]: 
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у/ ] = агсІ§[А 2 (со)/ 5, (со )} . (5) 

По рассчитанным значениям А ф ю, и <р ; - судят о 
спектральном портрете анализируемого сигнала 
/[{). Приведенные выше формулы (4) и (5) исполь- 
зованы в [1] для анализа спектрального состава 
сигнала в нелинейных системах автоматического 
регулирования. 

При переходе к решетчатым периодическим 
функциям формулы (3) принимают следующий вид: 

2 "м-і 

5|(ю,) = — I /( 4 ) 8 іп(ю. 4 ); 

М м к =° 

2 "м- 1 

5 2 (©,) = — I /(4)С08(Ю. 4 ); к 

1 Ѵ М *=о 

5(®.) = ^ 1 2 (ю.)+ 8 2 (®,), 

где 5(соу) - мгновенная спектральная плотность на 
определенной частоте опорного сигнала, /Ѵ ѵ - ко- 
личество точек на периоде анализируемого сигна- 
ла, попавших в моменты времени 4 при определен- 
ном М. При отсутствии в дискретизированом сиг- 
нале требуемых отсчетов из сигнала выбирают от- 
счеты, максимально приближенные к искомым 
моментам времени. 

Амплитудное значение у'-ой гармонической со- 
ставляющей рассчитывается по формуле: 

Лу= ^КО- 

Фазовый угол /-ой гармонической составляю- 
щей рассчитывается по формуле (5). Для иллю- 
страции работоспособности рассматриваемого 
способа спектрального анализа приведем расчет 
тестового примера с достаточно близкими друг к 
другу частотами, табл. 1, при шаге дискретизации 
А/=10 4 с и со мч = 2л50 рад/с, со пт = 2 л:5 1 рад/с, 
Дет=2 яО, 1 рад/с. Массив цифровых отсчетов сфор- 
мирован на периоде сигнала 7= 10 с. 


Таблица 1. Исходные данные тестового примера 


№ гармонической 

Частота 

Амплитудное 

Фазовый угол 

составляющей 

І=а>І2л, Гц 

значение Д,, В 

ср, градус 

і 

50,0 

200,0 

45,0 

2 

50,1 

200,5 

83,0 

3 

50,2 

200,3 

0,0 


Результаты вычислений тестового примера при 
различных значениях М представлены в табл. 2. 

По результатам экспериментов можно утвер- 
ждать, что рассматриваемый способ обеспечивает 
сокращение количества обрабатываемых отсчетов, 
и соответственно сокращает временные затраты на 
расчет спектрального портрета. 

Синхронный детектор 

Известный принцип синхронного детектирова- 
ния сводится к умножению анализируемого сигна- 
ла а(1) на опорный сигнал Ъ т (і) той или иной фор- 


мы, частоты и фазы (рис. 1). Среднее значение вы- 
ходного сигнала С(І) синхронного детектора четко 
зависит от амплитуды частотной составляющей 
сигнала а(і), причем при фазовой подстройке 
опорного сигнала удается выделить соответствую- 
щую составляющую сигнала практически без осла- 
бления [5]. 


Таблица 2. Результаты расчета тестового примера 


Способ 

Кол-во от- 
счетов, N 

Частота 
/2л, Гц 

Амплитудное 
значение А ш , В 

Фазовый 
угол ср , град. 

Класси- 

ческий 

МСП 

100000 

50,0 

200,0000 

45,0000 

50,1 

200,5000 

83,0000 

50,2 

200,3000 

0,0000 

при 
М= 4 

2000 

50,0 

200,0000 

45,0000 

2004 

50,1 

200,5000 

83,0000 

2008 

50,2 

200,3000 

0,0000 

при 
М= 5 

2500 

50,0 

200,0000 

45,0000 

2505 

50,1 

200,5000 

83,0000 

2510 

50,2 

200,3000 

0,0000 

при 
М= 6 

3000 

50,0 

200,000 

45,0000 

3006 

50,1 

200,5000 

83,0000 

3012 

50,2 

200,3000 

0,0000 

при 

М=1 

3500 

50,0 

200,0021 

44,7366 

3507 

50,1 

200,498 

82,9999 

3514 

50,2 

200,2954 

0,0000 

при 
М= 8 

4000 

50,0 

200,0000 

45,0000 

4008 

50,1 

200,5000 

83,0000 

4016 

50,2 

200,3000 

0,0000 



ЪоЛ О 

Рис. 1. Синхронный детектор 

В качестве опорного (зондирующего) сигнала 
рассмотрим решетчатый периодический сигнал 
вида Ъ оп (( к )=$іп(2л/ 0 „{ к + (//„„) =$іп(й)„„4+ уО, прини- 
мающий на своем периоде фиксированное количе- 
ство отличных от нуля значений: 

Ч я если (®«Л + ѵО=*/2; 

. , ѵ I -4. если (®„Л+ѵО = 3л /2; 

О если (, чЛ+ѵО^/2 , 
КЛ+ѵО = Зя/2. 

Опорный сигнал многократно сдвигают по фазе 
относительно анализируемого сигнала а(і) с шагом 
А у/ от 0 до 180°, причем частота опорного (зонди- 
рующего) сигнала со ж изменяется в пределах 
о) НАЧ <а> т <а> ПР щ. Находят точки (ю т 4+уО, в которых 
4,„(4) = Ля или _ Ч П Р И различных частотах опорно- 
го сигнала со 0 = 2л/ ж и фазах опорного сигнала у/ ж , 
перемножают значения анализируемого сигнала и 
опорного в точках Ь т ( 4), суммируют полученные 
значения на периоде анализируемого сигнала, рас- 
считывают постоянную составляющую, причем из 
найденных значений постоянной составляющей на 
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каждой частоте опорного сигнала выбирают макси- 
мальное значение С 0І и соответствующее ему значе- 
ние фазового угла опорного сигнала. Амплитудное 
значение выявленной гармонической составляю- 
щей рассчитывают по формуле: 

А ті =2(С ш /А т ). 

Среднее значение выходного сигнала С(1) будет 
максимально только при совпадении частоты и фа- 
зового угла опорного сигнала с частотой и фазовым 
углом реально присутствующей в анализируемом 
сигнале гармонической составляющей. Тогда, за- 
фиксировав положение опорного сигнала относи- 
тельно анализируемого сигнала при С срткс , можно 
определить фазу анализируемого сигнала. При от- 
сутствии в дискретизированом сигнале требуемых 
отсчетов, из сигнала выбирают отсчеты, максималь- 
но приближенные к искомым моментам времени. 

В табл. 3 представлены результаты расчета при- 
веденного выше тестового примера, табл. 1. Шаг 
перебора по частоте Д/ ш =0,1 Гц, по фазовому углу 
Ду/=Г. 


Таблица 3. Результаты расчета тестового примера 


Способ 

Кол-во от- 
счетов, N 

Частота 
Т^со/Ія, Гц 

Амплитудное 
значение А ш , В 

Фазовый 
угол <р , град. 

Класси- 

ческий 

ІО 6 

50,0 

200,0000 

45,0 

50,1 

200,5000 

83,0 

50,2 

200,3000 

0,0 

Предла- 

гаемый 

1000 

50,0 

200,0000 

45,0 

1002 

50,1 

200,4790 

83,0 

1004 

50,2 

200,2920 

0,0 


На рис. 2 представлен случай, когда анализиру- 
емый сигнал близок к синусоидальному, а опорный 
сигнал имеет единичную амплитуду при частоте, 
близкой к частоте основной гармоники анализиру- 
емого сигнала (Т 0 =Т С ). Тогда среднее значение вы- 
ходного сигнала при опорном сигнале Ь ш] (і) будет 
заметно больше, чем при опорном сигнале Ъ т1 

С . > С т, 

срі ср2 ’ 

что объясняется несовпадением фазы опорного 
сигнала у/ 2 с фазой анализируемого сигнала (в кон- 
кретном случае она равна нулю), тогда как для 
опорного сигнала Ь т] (і) 1//,=0. 

На рис. 3 рассмотрен случай, когда опорные 
сигналы отличаются друг от друга по частоте, при- 
чем Т т1 т с , т оп2 >т с . Легко видеть, что средние зна- 
чения выходного сигнала будут заметно отличаться 
друг от друга. Этот факт может быть использован 
при контроле частоты сигнала а(і) [8]. 

Отметим, что постоянная составляющая выход- 
ного сигнала С ср будет максимальна при выполне- 
нии двух условий, приведенных на рис. 2 и 3, если 
совпали: 

• частоты опорного сигнала и основной соста- 
вляющей анализируемого сигнала; 

• фазовые углы гармонической составляющей и 
опорного сигнала. 


По результатам вычислительных экспериментов 
видим, что применение решетчатых периодических 
функций дает значительное уменьшение числа об- 
рабатываемых отсчетов, то есть обеспечивается су- 
щественное уменьшение временных затрат. 

Дискретное преобразование Фурье 

Формулы дискретного преобразования Фурье 
хорошо известны [4]. При переходе к решетчатым 
функциям будем иметь: 

2 А ' м ~ 1 

Л = — X Л^) с °8( 2 ^^/л 2 ); 

Дм у=о 



А тк =уІ4+в I, 

'А* =агсІё(Л,/5,), 

где /|=1/ Г- основная частота сигнала, Т - период 
анализируемого сигнала, к - порядковый номер 
спектральной (гармонической) составляющей ана- 
лизируемого сигнала. 
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Рис. 3. Выходные сигналы синхронных детекторов при 
Т„„і = Т с и Т т1 > Т с 

Решетчатые периодические функции ^ й (/,) и 
принимают на своем периоде строго фикси- 
рованное количество значений в определенные мо- 
менты времени г, ( 1) при различных частотах гармо- 
нического сигнала о) к =ка> 1 =к2я/ ь кратных основ- 
ной частоте анализируемого сигнала/^ 1/ Г. 

Из анализируемого сигнала выбирают только те 
отсчеты, которые попадают в моменты времени 1 І 
при различных частотах со І =к(о 1 =к2л/ [ . Перемно- 
жение значений анализируемого сигнала и опор- 
ного осуществляется только в моменты времени ( к . 
Полученные значения суммируют на периоде Т 
анализируемого сигнала и рассчитывают значения 
А тк , со к и у/ к . При отсутствии в дискретизированом 
сигнале требуемых отсчетов, из сигнала выбирают 
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отсчеты, максимально приближенные к искомым 
моментам времени. 

Работоспособность расчетных формул и досто- 
верность получаемых результатов рассмотрим на 
тестовом примере, табл. 1. Результаты расчетов 
сведены в табл. 4. 

По приведенным в табл. 4 результатам расчетов 
видно, что предложенный способ расчета дискрет- 
ного преобразования Фурье с использованием ре- 
шетчатых периодических функций дает верные ре- 
зультаты расчета при различном М. 


Таблица 4. Результаты расчета тестового примера 


Способ 

Кол-во от- 
счетов, N 

Частота 
^со/2л, Гц 

Амплитудное 
значение Д,, В 

Фазовый 
угол (р , град. 

Класси- 

ческий 

МСП 

100000 

50,0 

200,0000 

45,0000 

50,1 

200,5000 

83,0000 

50,2 

200,3000 

0,0000 

при 
М= 4 

2000 

50,0 

200,0000 

45,0000 

2004 

50,1 

200,5000 

83,0000 

2008 

50,2 

200,3000 

0,0000 

при 
М= 5 

2500 

50,0 

200,0000 

45,0000 

2505 

50,1 

200,5000 

83,0000 

2510 

50,2 

200,3000 

0,0000 

при 
М= 6 

3000 

50,0 

200,0000 

45,0000 

3006 

50,1 

200,5000 

83,0000 

3012 

50,2 

200,3000 

0,0000 

при 
М= 1 

3500 

50,0 

200,0021 

44,7366 

3507 

50,1 

200,4980 

82,9999 

3514 

50,2 

200,2954 

0,0000 

при 
М= 8 

4000 

50,0 

200,0000 

45,0000 

4008 

50,1 

200,5000 

83,0000 

4016 

50,2 

200,3000 

0,0000 


В заключение обзора возможных путей приме- 
нения решетчатых периодических функций при 
расчете спектральных портретов узкополосных пе- 
риодических сигналов еще раз отметим существен- 
ное сокращение количества обрабатываемых от- 
счетов. Выбор оптимального значения М осущест- 
вляется экспериментальным путем при решении 
конкретных задач и зависит от доступного набора 
частот дискретизации анализируемого сигнала. 
Разработанные способы могут быть использованы 
при построении адаптивных систем управления с 
использованием принципа частотно-зависимой 
обратной связи с многочастотным узкополосным 
идентифицирующим воздействием. 

По результатам исследований были получены 
патенты [6-8], отдельные вопросы использования 
решетчатых периодических функций частично 
описаны в [3]. 
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